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Abstract. Using next-to-leading order perturbative QCD, we find that the charm 
contribution to the atmospheric lepton fluxes dominates over the conventional ones 
from 7r and K decays for energies higher than 10 5 GeV. We also discuss theoreti- 
cal uncertainties involved in the calculations and compare our results with previous 
evaluations. 



INTRODUCTION 

Conventional lepton fluxes from ir and K decays in the atmosphere are well known 
[1], but at energies of a few GeV charm- ant icharm pairs can be produced and the 
leptons coming from charm decays, called prompt leptons, start to contribute to 
the atmospheric lepton fluxes. Results from several experiments [2] show a muon 
excess relative to the conventional muon flux in the 10 TeV energy range that could 
be an indication of a charm contribution. Large volume water and ice neutrinos 
experiments [3] are designed to search for muon neutrinos from extragalactic sources 
for which atmospheric muons and muon neutrinos represent the main background. 
Therefore it is important to determine the atmospheric lepton flux, including charm 
contributions, with better accuracy. 

Atmospheric lepton fluxes from charm decays have been evaluated in the past 
using different models of charm production [4,5]. We present a new calculation 
[6] based on next-to-leading order (NLO) perturbative QCD. We also investigate 
the importance of the small-x behaviour of the parton distribution functions and 
study the uncertainties involved in the necessary extrapolation of cross sections 
and energy distributions beyond the experimentally measured regime. 



LEPTON FLUX CALCULATION 

The lepton fluxes can be calculated semi-analytically using approximate cascade 
equations which begin with the incident cosmic ray flux. Details of this solution 
method can be found, for example, in Refs. [7,6]. As in Ref. [4], we assume that 
the incident cosmic ray flux can be represented by protons and that at the top 



of the atmosphere (X=0) the proton flux (f) p (E,X) scales as 1.7 (E/GeV) 2 - 7 for 
£ < 5 • 10 6 GeV and scales as 174 (E/GeV)~ 3 for E > 5 • 10 6 GeV. 

Production of charmed hadrons j of energy E by proton-air interactions enters 
into the calculation via the "Z-moment", expressed in terms of an integral over 
Xe = E/E p , which is given by 

Z =2f f 1 dx E</> P (E/x E ) 1 da pA ^ c5 (E/x E ) 
Pj ~ Ij Jo x E <t> p {E) a pA (E) dx E ~ ~ U 

where fj is the fraction of charmed particles which emerge as hadron j and the 
factor of two is taking into account the fact that we implicitly sum over particles 
and antiparticles. The quantity g p a{E) is the total proton-air cross section [8]. 

The charm production cross section and energy distribution depends on the value 
of the charm mass (m c ), the factorization scale (M) and renormalization scale (//) as 
well as the choice of parton distribution functions. Comparing the NLO total charm 
cross section [9] with the experimental data from Ref. [10] led us to set m c = 1.3 
GeV. We use the CTEQ3 [11] and the MRSD- [12] parton distribution functions. 
We set our default factorization and renormalization scales to M — 2/j, — 2m c . We 
also evaluate the lepton flux using the CTEQ3 parton distribution functions with 
M = ji = m c . 

Because of the large computing time involved, we have used the NLO Monte 
Carlo program from Mangano, Nason and Ridolfi [13] to evaluate the appropriate 
rescaling factor K that, applied to the "LO" energy distribution, correctly accounts 
for the NLO corrections. Here "LO" means taking the leading order matrix element 
squared, but using the two- loop a s (/i 2 ) and the NLO parton distribution functions. 
The parametrization of K(E,x E ) is given by the following expression 

K(E, x E ) = 1-36 + 0.42 ln(ln(£/GeV)) (2) 
+ (3.40 + 18.7(£/GeV)- a43 - 0.079 ln(£/GeV)) • x E 5 . 

The i^-factor rescaled "LO" charm production distribution is used in Eq. (1) to 
obtain our lepton fluxes. 

In Fig. 1 we show our results for the prompt atmospheric muon flux scaled by 
E 3 for two parton distributions and factorization scale choices. The highest flux at 
E = 10 s GeV is with the MRSD- distribution and M = 2/x = 2m c (dashed). The 
CTEQ3 distributions with the same choice of scale are represented by the solid line, 
while the dot-dashed line shows the result when M — /j, — m c . For reference, we 
show the vertical conventional and prompt muon flux calculated and parameterized 
by Thunman et al. in Ref. [4]. The prompt flux is isotropic. The prompt lepton 
flux evaluated using perturbative QCD can be parameterized as 

log lo (£ 3 £ (£)/(GeV7cm 2 ssr)) = -A + Bx + C x 2 - Dx 3 , (3) 

where x = log 10 (£'/GeV). In Table I, we collect the constants for the MRSD- and 
CTEQ3 fluxes exhibited in Fig. 1. The prompt neutrino fluxes equal the prompt 
muon flux. 



TABLE 1. Parameters for the prompt muon plus antiunion 
flux appearing in Fig. 1, parameterized by Eq. (3). 



PDF Scales A B C D 

CTEQ3 M = fx = m c 5.37 0.0191 0.156 0.0153 

CTEQ3 M = 2fj,= 2m c 5.79 0.345 0.105 0.0127 

MRSD- M = 2[i = 2m c 5.91 0.290 0.143 0.0147 
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FIGURE 1. The prompt atmosheric muon flux scaled by E 3 versus muon energy for CTEQ3 
(solid) and MRSD- (dashed) with M = 2fi = 2m c . Also shown is the scaled muon flux using 
CTEQ3 with M = /j, = m c (dot-dashed) and the Thunman et al. parametrization [4] of the 
prompt muon flux and the vertical conventional muon flux (dotted). 




CONCLUSIONS 

We find that the perturbative charm contributions to lepton fluxes are signifi- 
cantly larger than the recent Thunman et al. calculation [4]. The prompt muon 
flux becomes larger than the conventional muon flux from pion and kaon decays at 
energies above ~ 10 5 GeV. We have evaluated the energy distribution at NLO in 
QCD, and we find that the NLO corrections give a correction of more than a factor 
of two which is weakly energy and xe dependent. 

The main uncertainty in the perturbative calculation of the prompt flux, given 
fixed charm mass, factorization scale and renormalization scale, is the small-x be- 
havior of the parton distribution functions. The spread in predictions at E — 10 8 
GeV is indicative of this uncertainty. 

We conclude that the prompt muon flux calculated in the context of perturbative 
QCD cannot explain the observed excess of muons in the TeV region [2], indepen- 
dent of the theoretical uncertainties associated with small parton x. Measurements 
of the atmospheric flux in the 100 TeV range would help pin down the charm cross 
section at energies above those currently accessible using accelerators and would 
provide valuable information about the small- x behavior of the gluon distribution 



function. 
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